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Abstract: In this paper we investigate the localization and mass spectra of bulk matter 
fields on a Gergory-Rubakov-Sibiryakov-inspired braneworld. In this braneworld model, 
there are one thick brane located at the origin of the extra dimension and two thin branes 
at two sides. For spin 1/2 fermions coupled with the background scalar <f> via r]fy(f> p ty 
with p a positive odd integer, the zero mode of left-hand fermions can be localized on the 
thick brane for finite distance of the two thin branes, and there exist some massive bound 
modes and resonance modes. The resonances correspond to the quasi-localized massive 
fermions. For free massless spin scalars, the zero mode can not be localized on the thick 
brane when the two thin branes are located finitely. While for a massive scalar <I> coupled 
with itself and the background scalar field 4>, in order to get a localized zero mode on the 
thick brane, a fine-tuning relation should be introduced. Some massive bound modes and 
resonances also will appear. For spin 1 vectors, there is no bound KK mode because the 
effective potential felt by vectors vanishes outside the two thin branes. We also investigate 
the physics when the distance of the two thin branes tends to infinity. 
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1. Introduction 

Recently, the braneworld theory has attracted more and more attention [jl], 0, ||, ||, [5l ||] . In 
this theory, our matter fields are localized on the 3-dimensional brane, which is embedded 
in a higher dimensional space-time, while gravity is free to propagate in the whole space- 
time. The most important application of this theory is that it can solve some particle 
problems, such as mass hierarchy problem and cosmological constant problem. 

Two famous brane models brought forward in early times are the Arkani-Hamed- 
Dimopoulos-Dvali (ADD) large extra dimension brane model Q and the Randall- Sundrum 
(RS) warped brane model [||. The extra dimensions of the ADD model are flat and 
compact, and the ADD model predicts the deviations from the 4D Newton law at sub- 
millimeter distances. While the RS model has an infinite extra dimension using a non- 
factorizable "warped" geometry. And in the RS model, there is a massless bound state 
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graviton which should be interpreted as the usual 4D graviton. However, they are both 
thin branes, which have no thickness. 

On the base of the RS model, another model called the Gergory-Rubakov-Sibiryakov 
(GRS) model was proposed 0, in order to give some modifications of gravity at ultralarge 
scales. There are three thin branes in this model, which are located at y = 0, y = yo 
and y = —yo (hereafter y denotes the fifth coordinate), respectively. In the Ref. ||], the 
GRS-inspired brane model was set up, in which the central thin brane is replaced by a 



thick brane which is realized by a scalar field, just as the usual thick branes |9j, |10j, [llj, |12 

n n m m, mm. 

There are also other braneworld models, but for all models, an important problem is 



the localization of various matter fields on the brane [20, 21, 22, 23, 24]. It has been known 
that massless scalar fields and graviton can be localized on branes of different types 0, . 
But spin 1 Abelian vector fields can not be localized on the RS brane in five dimensions, 
however, it can be localized on the RS brane in some higher-dimensional cases |26|], or on 
the thick de Sitter brane and the Weyl thick brane p7|, |28|] . 



The localization of fermions is more interesting [29, 3C, 31, 32, |33|] and has been ex- 



tensively analyzed in many papers, such as with gauge field background [34, 35], with 
supergravity background [36], with vortex background [37, 38, 39, 40 [. It was also found 
that there may exist a single bound state and a continuous gapless spectrum of massive 
fermion Kaluza-Klein (KK) states gl|, g|, [|| (£| ||, g|], or finite discrete KK states 
(mass gaps) and a continuous spectrum starting at a positive m 2 [28, 47, 48] with the 



scalar- fermion coupling. Some resonances also exist [49, 50 1, whose life-time depends on 
the parameters of the scalar potential. 

In this paper we investigate the localization and mass spectra of bulk matter fields on 
the GRS-inspired braneworld. Our paper is organized as follows: In Sec. ||, we first give a 
brief review of the GRS-inspired braneworld. Then, in Sec. ^, we study the localization and 
mass spectra of spin 1/2 fermions, spin scalars and spin 1 vectors on the GRS-inspired 
braneworld by presenting the shapes of the potentials of the corresponding Schrodinger 
equations. Finally, in Sec. ||, we make a brief discussion and a conclusion. 



2. Review of the GRS-inspired model 

An interesting model which has the feature of predicting modifications of gravity at ul- 
tralarge scales is the one proposed by Gregory, Rubakov and Sibiryakov (GRS). In the 
usual GRS model, considering the fifth dimension y to be infinite and non-compact, there 
are three thin branes, one is at y = with positive tension a, the other two are at y = yo 
and y = —yo, respectively, both with the tension —a/2. Note that the total charge in this 
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braneworld is zero. The cosmological constant A between the positive and the negative 
tension branes is assumed to be negative, namely, the space for \y\ < yo is AdS$, which 
is the maximally symmetric space. Beyond |yo| the space is the Minkowski space. The 
system is assumed to have orbifold symmetry. So we can deal with the domain y > only, 
because for negative values of y the system is just a mirror reflection. 

If the central thin brane is replaced by a thick brane realized by a suitable configuration 
of a bulk scalar field which is called the GRS-inspired brane Q, the action of this 
model is 

S= I d':v^ -G 



2M 3 R - ^(V0) 2 - V(<j>) - A 



J d^xy/^ga^), (2.1) 



where M is the basic Planck energy scale in the 5-dimensional space-time, Gmn denotes 
the 5-dimensional metric and g^ v the 4-dimensional one. The first part of the action 
contains the gravitation and the scalar field, while the second part is the contribution of a 
thin brane. The background scalar <p is considered to be only the function of the extra 
dimensional coordinate for easier, i.e., <j) = 4>(y), and the brane tension depends on the 
scalar field only, i.e., a = a ((f)), because the branes are treated as the cross-sections of the 
bulk. 

In general, the geometry of a 4-dimensional brane with maximal symmetry can be 
supposed as M4, dS^ or AdS^. Here, we consider the simplest case, the M4 brane. Hence, 
the metric of the background space-time is assumed as |K 



ds 2 = e 2A{ - y \^ v dx^dx v + dy 2 , (2.2) 
where e 2A ( y "> is the warp factor. For convenience, the metric (2.2) can be rewritten as 



ds 2 = e 2A(z) (r/^dx^dx" + dz 2 ) , (2.3) 

where the two coordinate systems are connected by dy = e A ^dz, so the point y = is 
mapped to z = 0, while yo to z$. And the extra dimension z is also infinite. 
Then the equations of motion generated from the action ( |2.1| ) are 

Rmn — -^GmnR 



2 4M 3 



Va/^Vtv^ - GW(^(V0) 2 + V (</>)) 



-Gmn A - ct8(z - z )e~ A ^ 9llu S M S^ 



(2.4) 



r2, dV 



W+^ = 0, (2.5) 
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which can be reduced to the following coupled nonlinear differential equations with the 
ansatz (|2.3|): 

<p' 2 = 12M 3 (A' 2 - A") - a5(z - z )e" A(2o) , (2.6) 
V(<f>) = 6M 3 {-3A' 2 - A")e- 2A - A - °-b{z - zo)e" A(2o) , (2.7) 

^ = + <P")e~ 2A ~ ^6(z - z )e-^\ (2.8) 

a<p cup 

where the prime stands for the derivative with respect to z. Note that, here we only consider 
the range of z > 0, so terms containing 5{z + zq) do not appear in above equations. 
The background solution for this GRS-inspired domain wall is H 1 

• 4 = i ln rnw' (2 - 9) 

<p = v / 12M 3 / 2 arctan(A;z), (2.10) 
.■■^f ^'-V i. (2.11) 



12M 3 / 2 - 



for \z\ < zq, and 

A=-hx r — ttt, (2.12) 

(p = \/l2M 3/2 arctan(A;zo), (2.13) 
V = -A, (2.14) 

for \z\ > zq, where k and b are both constants and k plays the role of the central wall's 
thickness. With the normalization condition e 2A ^ = 1 for the warp factor, the parameter 
b can be set to 1. And with the solution of the warp factor above we can calculate the 
exact relation between the two coordinates: 

rz 

y = J e A(5) dz = -arcsinh(fcz), (2.15) 

for — zq < z < zq, and 

y = f e A ^dz = c + ? =L= z, (2.16) 
Jo V 1 + k 2 zfi 

for \z\ > zq, where c is a constant as c = iarcsinh^Zo) H — . z ° =. . So it is clear the extra 
coordinate z is also infinite. 



lr The potential 1/ in (2.11) for z > zq is a little different from that in Ref. fl8|, where V 



L- 2 A/3 



-A. 
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It can be seen that V(cfi) may be discontinuous at zq, unless we set zq = I /{2k). 

v _ 3fc 2 [5(fcz) 2 -2] ■ 
Vl ~ 4[l+(fcz) 2 ] 2 



However, the localizing potential of gravity for \z\ < zq, i.e., V\ 



can not 



reach zero if we set zq = ^ because the zero point of V/, Z\ = yfp is larger than 
So, considering that Vi = for \z\ > zq (this is because we have a Minkowski space for 
\z\ > zq), we would get a localizing potential with no barrier if we perform the patch. The 
detailed discussion can be found in Ref [p. 

The scalar configuration is in fact a kink. It is clear that beyond zq the warp factor 
and the scalar field are both constants. So the expression of the energy density p is 



12fc 2 (l-2fc 2 z 2 ) 
l+k 2 z 2 



\Z\ < Z 

f(l + k 2 zlfl 2 [5(z-z ) + d(z + z )], \z\ >z 



(2.17) 



The tension of the thin branes a satisfies a 



l2M 3 k : 



which is calculated from the 



junction condition at z = ±zq |8|. We can see that the energy density vanishes at \z\ > zq, 
which is because the space beyond zq is Minkowski. For convenience, we set M = 1 in 
what follows. The shapes of the warp factor e 2A ( z \ the kink 4>(z), and the energy density 
p(z) are shown in Fig. [l] and Fig. ^. 




Figure 1: (Color online) The shapes of the warp factor e 2A ^ and the kink <f>(z). The parameters 
are set to k = 1, Zq = 0.5 for red thick line (which is the case 2fczo = 1), zo = 3 for green dashing 
line, and zq — 6 for blue thin line. 



3. Localization and mass spectrum of various matter fields on the brane 

In this section, we will investigate the localization and mass spectra of various bulk matter 
fields such as spin 1/2 fermions, spin scalars and spin 1 vectors on the GRS-inspired 
braneworld by presenting the potentials of the corresponding Schrodinger equations. We 
treat the bulk matter fields considered below as small perturbations around the background 
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Figure 2: (Color online) The shapes of the energy density p(z). The parameters are set to zq = 5, 
k = 0.1 for the left figure (which is the case 2kzo = 1), and k — 3 for the right figure. 



1 53, 54], namely, we neglect the back-reaction of bulk matter fields on the background 
geometry. So the bulk fields make little contribution to the bulk energy, and the background 
solutions given in previous section are still valid. 

3.1 Spin 1/2 fermion fields 

Firstly, we will investigate the localization and mass spectra of the spin 1/2 fermions on 
the GRS-inspired braneworld by means of the gravitational interaction and scalar-fermion 
coupling. 

Using our usual set-up, we consider a massless spin 1/2 fermion coupled with gravity 
and the background scalar in 5-dimensional space, of which the Dirac action is 

S 1/2 = Jd 5 xV^G(^T M (d M (3.1) 

where T M = (e~ A j^ 1 ,e~ A 'j 5 ), uj^ = 5(^.4)7^75, 0J5 = 0, and 77 is the coupling constant 
between the spinor ^ and the scalar field eft with F(cft) the type of the coupling. 7^ and 7 s 
are the usual flat gamma matrices in the Dirac representation. The 5-dimensional Dirac 
equation is read as 

{Vfy + 7 5 (d z + 2d z A) - r, e A F(0)} V = 0, (3.2) 

where 7^6^ is the 4-dimensional Dirac operator. 

Now we would like to decompose the 5-dimensional spinor ^ as 

(Z, Z) = e~ 2A ULn(x)f L n(z) + ^Rn(x)f Rn (z)\ , (3.3) 
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where tpLn = —J 5 ipLn and ipR n = ^ 5 i^R n . Substituting the above expression of ^(x, z) into 
the equation of motion (3.2), and assuming that ipLn,Rn(x) satisfy the 4-dimensional mas- 



sive Dirac equations j fl d^L nt R n (x) = rn n tpR nt Ln(x), one can obtain the coupled equations 
of the fermion KK modes fin(z) and fR n (z): 

[d z + r, e A F X<f>)] hn(z) = m n f Rn (z), (3.4a) 

[d z - rj e A F(ct>)] f Rn (z) = -m n f Ln (z). (3.4b) 

Then from the above coupled equations, we get the Schrodinger-like equations for the KK 
modes of the left and right chiral fermions: 

( - d 2 z + V L {z))f Ln = m 2 J Ln , (3.5a) 

(- d 2 z + V R (z))f Rn = m 2 J Rn , (3.5b) 

where the effective potentials are 

V L (z) = {ve A F(</>)) 2 -d z {r ] e A F{ct>)), (3.6a) 

Vr(z) = V L (z)\ v ^ n . (3.6b) 

Moreover we need the following orthonormality conditions for fm and fn n to get the 
standard 4-dimensional action for massive chiral fermions: 



fLm{z)fLn{z)dz = 5 mn , 

) 

fRm(z)fRn(z)dz = 5 mn , (3.7) 

) 

fLm(z)f Rn (z)dz = 0. 

) 

It can be seen that the potentials Vl(z) and Vr(z) depend on three factors, i.e., the 
scalar-fermion coupling constant rj, the warp factor and the type of scalar-fermion 

coupling F ((/)). In order to localize left- or right-hand fermions, the effective potential Vl(z) 
or Vr(z) should have a minimum value at the location of the thick brane. Furthermore, 
as we demand that Vl,r(z) are invariant under Z2 reflection symmetry z — > —z, F(cf>(z)) 
should be an odd function of z. Since the kink configuration of the scalar <f>(z) is ^2-odd, 
F(cj)) should be an odd function of (p. Thus we have F((f>(0)) = and Vl(0) = — Vr(0) = 
— r]d z F ((f)(0)), which would result in that only one of the massless left and right chiral 
fermions can be localized on the brane. 

From the solutions of the metric and the scalar field given in previous section, we can 
see that, for any type of scalar-fermion coupling, when \z\ > zq the value of the potentials 
is a constant, i.e., Vl : r(±oo) = Vl,r(zq). So there may exist discrete bound KK modes 
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for < Vl,r(zo) and continuous ones for > Vl,r(zo). Here, for simplicity, we only 
consider the simplest Yukawa coupling F(<j>) = <j) and the generalized one F(<p) = (jf with 
p an odd positive integer. 

For the simplest Yukawa coupling F{4>) = (ft, the potentials take the following form 



. i+fc^ + (1+k 2 z 2 )3 /2 L^arctan(fczj - Ij , \z\ < z 



Vr(z) = V L \ V ^ V , (3.8b) 
and for the generalized Yukawa coupling with p > 3, the potentials are 

J 12Pn 1 +° fc t ^ P(fcZ) + (i+k^ z 2)3/2 ^ctan p - 1 (kz)[kzaictan(kz) -p], \z\ < zq 

[ l+fc% 2 ' N ^ z o 

Vr(*) = V L \ V ^ V . (3.9b) 

The values of the potentials for left and right chiral fermions at z = are 

V L (0) = -V R (0) = -2V3kr), for p = 1, (3.10) 
Vr(0) = -Vk(0) = 0, for p > 3, (3.11) 

and V^^ioo) = Vl,_r(-< ; o) at y — > ±00, so it is clear the potentials are Poschl- Teller-like 
potentials. For positive coupling constant, the minimum value of the potentials for left- 
hand fermions is negative and for right-hand fermions positive, so only left-hand fermions 
have one zero mode. While for massive KK modes, the mass spectra are the same for both 
chiral fermions, which can be seen from ( p. 4a ) ( 3.4b| ). Because Vl,r(z > zq) = Vl,r{zq) is 
a constant, so the value of the potentials at z = zq is very important for the number of 
the massive bound KK modes. The value of Vl r{zq) is determined by three parameters: 
zq, r] and p. It can be seen that that Vl,r(zq) increases with r\ and p. And there is a 
special point of zq, for which Vl,r(zq) takes its maximum (see Fig. ||). The shapes of the 
potentials for different parameters are plotted in Fig. ||, Fig. ^| and Fig. || with Fig. |3| and 
Fig. [| for the simplest Yukawa coupling, and Fig. |5| for the generalized one. 

3.1.1 The zero mode 

The zero mode for the left-hand fermions can be solved from ( |3.4a ) by setting uiq = 0: 

f L0 (z) oc exp (-r)J* dze A &0(z^j . (3.12) 

In order to check the normalization condition ( |3.7|) for the zero mode ( p. 12 ), we need to 
check whether the integral 

J fl (z)dz oc J exp (-2ri J dze A ^<j){z^j dz (3.13) 
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Figure 3: (Color online) The shapes of potentials Vl (the left figure) and Vr (the right figure) 
with F((f>) — cj> for different coupling constant 77. The parameters are set to k = 1, zq = 5, 77 = 3 for 
red thick lines, r/ — 5 for green dashing lines, and rj = 7 for blue thin lines. 
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Figure 4: (Color online) The shapes of potentials Vl (the left figure) and Vr (the right figure) 
with F((j)) = (f> for different zo. The parameters are set to fc = 1, 77 = 3, zo = 0.7 for red lines, zq = 2 
for green dashing lines, and zq = 3 for blue thin lines. 



is finite. For the integral J z dze A (f>, we only need to consider the part of \z\ > zq, because 
the integral is finite for the integrand range — zq<z<zq. Since 



2e y 



12 



arctan(/czo) = c for \z\ > zq, 



V 1 + (^o) 2 

where c is a positive constant for positive k, the integrand f1 is 

/l (z)oce- 2 ^ for \z\>zo, 



(3.14) 



which indicates that for a finite zq the integral ( |3.13| ) is finite and the zero mode of left- 
hand fermions can be localized on the brane without any other condition. However, for an 
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Vl.Vr 




(a) p = 3 (b) p = 5 



Figure 5: (Color online) The shapes of potentials Vl (the red thick dashing lines) and Vr (the 
blue thin lines) with F((f>) = 4> p for different p. The parameters are set to k = f , zo = 3, ry — 0.02, 
p = 3 for the left figure, p — 5 for the right figure. 

infinite Zq, the integral J 2 dze A (f) tends to: 

/z ^/3tt 
dze A (j) — > ^ In z, for y — > ±oo, (3.15) 

then the integrand /| becomes 

/goW«z— r^, (3.16) 

so only when ri > 7L , the zero mode for the left chiral fermion can be localized on the 
brane, which is corresponding to that in Ref [plf| . 

3.1.2 The massive bound modes 

About the massive KK modes for both left-hand and right-hand fermions, as discussed 
before, the only difference between the mass spectra of left-hand and right-hand fermions 
is that there exists one zero mode for left-hand fermions but does not for right-hand ones 
for any positive coupling constant rj. In the following discussions we will see the impact of 
different parameters on the KK modes through their impact on the potentials. 

From Fig. ^, we can see that with the increase of 77, the depth of the potential well for 
both chiral fermions increases, which means that more massive bound KK modes would 
appear. Using the numerical method, we get the mass spectrum of the KK modes for both 
chiral fermions with k = 1, zq = 1 and different 77. For 77 = 1, there is only the zero mode 
for left-hand fermions. For 77 = 3, there are two massive bound KK modes for both chiral 



- 10 - 



fermions. For 77 = 5, there are three ones. The mass spectra are 



m \ n = {0, 18.19, 30.64} U [33.01, 00) for rj = 3, 
m 2 Rn = { 18.19,30.64} U [33.01, 00) for r) = 3, 

m\ n = {0, 32.09, 58.93, 79.97} U [92.53, 00) for 77 = 5, 
m 2 Rn = { 32.09,58.93,79.97} U [92.53, 00) for rj = 5. 



(3.17) 



(3.18) 



Next let us discuss the effect of the type of the scalar-fermion coupling, i.e., the effect 
of p. From Fig. ||, it can be seen that with the increase of p, the depth of the potential well 
for both chiral fermions increases sharply, which means that there would be more bound 
KK modes. Here we calculate the mass spectrum for p = 3, k = 1, rj = 1 and zq = 1: 

m \ n = {0,10.88,33.47,58.58,86.42, 116.00, 146.60, 177.44} U [202.79, 00), 

(3.19) 

m Rn = { 10.88, 33.47, 58.58, 86.42, 116.00, 146.60, 177.44} U [202.79, 00). 

Compared with the simplest Yukawa coupling (the case with p = 1), it is clear that when 
the coupling index p increases, the potentials could trap more bound KK modes. 

Finally, we come to the parameter zq, i.e., the location of the thin brane. Since the 
mass of a bound KK mode must satisfy the relation m n < \/Vl,r(zq), the effect of the 
location of the thin brane on the number of bound KK modes would be very important. We 
plot the shape of Vl(zq) as the function of zq in Fig. |6|. It is clear that, for the generalized 
Yukawa coupling, there is a maximal value of the potential. Using the first derivative of 
Vl{zq) about zq, we obtain that when the formula 

p = kzo arctan(A;2!o) (3.20) 

is satisfied, the value of Vl(zq) is maximum, which means the number of the bound 
KK modes is the most. For the simplest Yukawa coupling we found the formula is 
1 = kzQ arctan(/czo), meaning when kzo = 1.16, the value of Vl{zq) is maximum. Here 
we present the mass spectra with different zq for the simplest Yukawa coupling: 

m \ n = {0, 18.19, 30.69} U [34.00, 00) for z = 1.16, 

m\ n = {0, 16.90} U [18.57, 00) for z = 0.5, (3.21) 

ml n = {0 }U [2.31,oo) for z = W, 

where we set k = 1, rj = 3. It is clearly that only when zq = 1.16, the number of the bound 
KK modes is the most. If z$ — > 00, there will be only a zero mode for left-hand fermions, 
as V LtR (z ) -)• 0. 
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Figure 6: (Color online) The shape of the potential of left-hand fermion KK modes Vl(zq) (as the 
function of zo). The parameters are set to k = = 3,p = 1. 

3.1.3 The resonance modes 

Next, we also investigate quasi-localized KK modes of fermions, i.e., the fermion resonances. 
From the shapes of the potentials in Fig. |3| and Fig. we can see that when the position 
of the thin brane zq is not too small and the coupling constant r] is large, the potentials 
Vlr{z) would have a potential barrier at each side of the location of the thick brane. Those 
KK modes with m 2 < Vl,r(zq) are bounded and discrete, while those with m 2 > Vl,r(zq) 
are non-bounded and continuous. The low-energy continuous modes will generically have 
very small amplitudes at the position of the thick brane because they must tunnel through 
the potential barrier. Therefore, although low-energy fermions with m 2 > Vl,r(zq) would 
exist in the theory, their wave functions will be strongly suppressed on the thick brane, 
which leads to a weak coupling with the zero mode and massive bound modes. However, 
for higher potential barrier, the KK modes with certain discrete mass would resonate with 
the potential and will have larger probability of being found on the thick brane. 

From ( p.4a|) and fl3.4b|) we could conclude that the masses of the resonances for both 
chiral fermions are the same but their parity is reverse, which was discussed in detail in Ref 



1 55 J - Here we also use the method in that paper to study the resonances. First, in order to 
define the parity of the resonances, we give the initial condition, i.e., /(0) = h, /'(0) = for 
the even parity resonances, and /(0) = 0,/'(0) = h for the odd parity resonances, where 
h is a constant. Then we use the following expression to define the relative probability of 
the KK mode within a narrow range — Zj, < z < Zf, around the thick brane location []5l| 

f Zb \fL R {z)\ 2 dz 

w = ir iiiW (3 ' 22) 

where we choose z max = lOzj,. The large relative probabilities would indicate the existence 
of resonances. Here we give two examples, for p = k = 1, Zq = 5, rj = 3 and r] = 5. 
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We found two resonances for both chiral fermions for p = k = 1,Zq = 5, rj = 3, and 
the corresponding masses are given by m 2 = 18.1855 and m 2 = 30.60. The probabilities 
and the resonance wave functions are shown in Fig. and Fig. |8[ For p = k = 1,Zq = 
5,r/ = 5, there are four resonances for both chiral fermions, of which the masses are 
m 2 = 32.092, 58.93, 79.99, 93.7, respectively. Comparing with the massive bound KK 
modes with masses m 2 = 32.092, 58.93, 80.00, 93.7, where the parameters are set p = k = 
1, zq = 1.16, r] = 5, we can see that, when the distance of the two thin branes increases, the 
discrete massive bound modes would disappear, and the resonances are the remnants of 
the bound modes. Note that any particle with m 2 > Vl,r(zo) produced on the brane can 
not correspond exactly to a single resonant KK mode because it has a wave function truly 



localized on the brane [56]. It is a wave packet composed of the continuum modes with 
a Fourier spectrum peaked around one of the resonances. The quasi- localized resonances 



and fermion tunneling rate were investigated in [57] 



1.0 - 



0.4 - 



0.2 - 



Figure 7: (Color online) The probability Pi(m), Pfi(m) for finding the resonances of left and right 
chiral fermions around the thick brane location. The parameters are set to p = k = 1, r] = 3, Zq = 5. 
The red thick lines are for the even resonance, and the blue dashing lines are for the odd resonance. 



3.2 Spin scalar fields 

In this subsection, let us investigate the localization and mass spectra of spin scalars on 
the GRS-inspired brane. The action of a real scalar field <1> in 5-dimensional space-time is 

S = J d 5 x^G^~ ^G MN d M <S>d N <f> - F(<M)J , (3-23) 

where the first part describes the coupling of the scalar field ^ with gravity, while the 
second part V r ( < I>, (ft) is the coupling of the scalar field $ with itself and the background 
scalar field cj>, so it should include <I> 2 , <I> 3 , $ 4 and (4>&) 2 terms, but considering the 
discrete symmetry we have chosen, $,$ 3 terms can be eliminated. Note that, since the 
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Figure 8: (Color online) The shapes of the resonances /l and Jr. The parameters are set to 
p = k = 1,7} = 3,Zq = 5. The left pictures correspond to the resonances with m 2 = 18.1855, and 
the right pictures correspond to the resonances with m 2 = 30.60. 

scalar is introduced as a perturbation field, high order terms such as 3> 6 can be neglected. 
We can obtain the equation of motion using the conformally flat metric (|2.3|) 

d M (V z GG MN d N $) + V^G U(<f>)$ = 0, (3.24) 

where U is independent of and is denned as 9V ^'^ = U<f> + 0(<I> 2 ) [EjjL namely, 
we mainly consider the $ 2 and (</><3?) 2 terms in V(&,(j)). By decomposing <I>(x,,z) = 
Yin ^n(s ; )Xn(2 : )e~ 3j4 / 2 , we find the fifth part of the scalar field Xn{ z ) satisfies the following 
Schrodinger-like equation 

[-8 2 z + V (z)] xn(z) = m 2 n xn(z), (3.25) 
where we have assumed that the 4-dimensional massive Klein-Gordon equation 

(-^d^^tjg^dv) - m 2 ) <f> n (x) = (3.26) 
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is satisfied with g^ u = r]^ u the 4-dimensional metric on the brane, and the effective potential 
is 

Vo{z) = \^ + \^ 2 + e 2A U{4>). (3.27) 



If we take the scalar potential as the form [p6|, 55] 



0) = ( -M 2 - \u* ) & + 0(<I> 4 ), (3.28) 



^4 r 2 

where A and u are the coupling constants, then U(4>) can be read as 

U{4>) = \^4> 2 ~ u 2 - (3.29) 

Furthermore, we need the following orthonormality condition to obtain the 4-dimensional 
effective action: 

/DO 
dz Xm(z)Xn(z) = S mn . (3.30) 
-oo 

3.2.1 Free massless scalars 

We first consider the case without coupling, i.e., 5-dimensional free massless scalars with 
V($, 4>) = 0. The effective potential ( ^27|) is 

(- 3fc 2 [5(fcz) 2 -2] I I 

Fo( ^ = 4ii+(fo)T ' |Z| < Z ° . (3.31) 
I 0, \z\ > z 

We note that, if one replaces k 2 with 3A, the above potential between the two thin branes 



is just the one given in Ref. [44|, where the corresponding potential for scalar KK modes 
is Vq(z) = 9 4(i+3Xz 2 ~P • Hence, if we remove the two thin branes to infinity, we will recover 



the same result as in Ref. [44]. The zero mode can be solved as 



/ (1 + (MT 3/4 , \z\<zo 

Xo(z) K 1 , (i ^2^-3/4 , , . • ( 3 - 32 ) 
I 1 + {kzoY) ,\z\> z 



Clearly, the normalization condition (3.30) for the scalar zero mode is not satisfied for finite 
zq, which indicates that the 4-dimensional massless scalar can not be localized between the 
two thin branes with a finite distance. In fact, this result can be supported by the shape 
of the potential ( p. 31 ): it vanishes beyond the two thin branes, which indicates that there 
do not exist any bound KK modes including the zero mode. 

However, there are two limit cases we are interesting in: zq — > oo and zq <C 1/k. For 
the first case, the two thin branes are infinitely far away from the central thick brane, and 
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the effective potential ( 3,31 ) is in fact a volcano type potential. The zero mode Xo( z ) is 
normalizable and is turned out to be 



Xo{z) 



(3.33) 



V2(i + (k z yf 4 ' 

Namely, when we remove the two thin branes to infinity, we can localize 4-dimensional 
massless scalar on the thick brane. In addition to this massless mode, the potential ( 3.31| ) 
with zq — > oo suggests that there exists a continuum gapless spectrum of KK modes with 
positive m 2 > 0. 

Now we come to the second limit: zq <C 1/k, i.e., the distance of the two thin branes 
is much less than the thickness of the thick brane and the thick brane between the two 
thin branes has a constant density along the extra dimension. To the lowest order of 
approximation of kzQ, the effective potential (3.31) is the known finite square potential 
well: 

( 3 ,k 2 , \z\ < Z 
\z\ > Zq 



Vq{z) 



2' 

0, 



(3.34) 



Then we would get one and only one bound state with negative eigenvalue nig, which is 
decided by the relation: 



ml 



9 



-k 2 + m\ ) tan 2 



k 4 ZQ. 



(3.35) 



Hence, it seems that the system contains a tachyonic state in the limit of kzQ <C 1. However, 
if one use a more accurate potential (to the next order of kzo), i.e., 



Vo(z) 



\k 2 + ^k 2 Z 2 , \z\ < Zq 



0. 



\z\ > Zq 



(3.36) 



the result will be corrected. From the potential (3.36), we know that the corresponding 
solution of the bound KK modes is 



Xn(Z) 



Vbz 



r-jr: 2 — 'Hn^^z), \z\ < z 

~ m " z , Z > Zq 



c 2 e 



(3.37) 



c 3 e" 



z < -zo 



where H n (z) is the Hermite polynomial and n = (2m 2 + 3k 2 — 3\^3k 2 ) / \6\^3k 2 ) is an 
nonnegative integer. So the solution of the spectrum m 2 is 

9 3 



m r . 



-k 2 (2V3n + V3-l), 



(3.38) 
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from which we see that the eigenvalue is positive for any nonnegative n, and hence the 
potential ( p. 36 ) has no any bound state with negative eigenvalue, namely there do not 



exist any unstable KK modes . 



As for massive KK modes, we can use the approximate potential (3.34) to get the 
resonant spectrum of the scalar KK modes: 



n 2 7r 2 3k 2 n 2 7r 2 
= 1?— 2" * 1^" fr^ 1 ' 2 ' 3 '-)- (3-39) 



It is the same as the eigenvalue spectrum of the corresponding infinite square potential 
well. 

3.2.2 Scalars coupled with itself and the domain-wall-generating scalar <j) 

Next, we consider the case where the scalar <1> couples with itself and the domain- wall- 
generating scalar <f> with the scalar potential (|3.28| ). For the GRS-inspired brane solution, 
the effective potential (p.27|) can be expressed as 



2lk ' z2 + tto (6A arctan 2 ^) - \k 2 - u 2 ) , \z\ < z 



1+k 2 z i (oAaretan \kZq) — u ) , \z\ > zq 

We note that the value of the potential at the location of the thick brane is Vo(0) = 
— \k 2 — u 2 , which is always negative for nonvanishing k. While for \z\ > zq, the value 
of the potential is a constant Vo(zq). So when the relation 6A < u 2 arctan _2 (/czo) is 
satisfied, there will bound KK modes with negative eigenvalues m 2 because Vb(zo) < 0. 
Thus in order to get stable solutions, we need to introduce a fine tuning condition 6A = 
(u 2 + v 2 ) arctan~ 2 (A:zo) with v some constant to be decided by the parameters u, k, and zq, 
for which Vo(^o) ^ ' t;2 /(l + k 2 z 2 ) and the potentials become Poschl- Teller-like potentials. 
So with the fine tuning condition, there will exist discrete bound KK modes, resonances 
modes and continues non-bound KK modes with nonnegative eigenvalues and we will get 
a localized massless scalar on the thick brane. Moreover, in order to localize the massless 
or massive scalars on the thick brane, the coupling constant A should be large enough. 

The potential is impacted by three parameters, i.e., the coupling constant A, the mass 
parameter u, and the position of the thin brane zq. The shapes of the potential are shown 
in Fig. m for different parameters. 

According to the above analysis, we know that there exists a localized massless scalar 
on the thick brane if and only if a fine tuning relation is satisfied. This can be verified by 
numerical calculations. As a example, we found that when k = 1, Zq = 3, A = 23, u = 3.1215 
(the corresponding value of v is 14.3371), there exists a massless bound KK mode localized 
on the thick brane. The wave function of the zero mode is shown in Fig. |i~C|. 
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Figure 9: (Color online) The shapes of the effective potential Vo(z) of the scalar KK modes. In 
the left figure, the parameters are set to k — 1, zo = 3, A = 23, u — for red thick lines, u — 3 for 
green dashing lines, and u — 5 for blue thin lines. In the middle figure, the parameters are set to 
k = 1, Zo = 3, u = 1, A = 23 for red thick lines, A = 43 for green dashing lines, and A = 53 for blue 
thin lines. In the right figure, the parameters are set to k = 1, A = 43, u = 1, Zq = 0.5 for red thick 
lines, zq = 1.16 for green dashing lines, and zq = 3 for blue thin lines. 



Xo 




-2-1 12 



Figure 10: (Color online) The shape of the zero mode xo f° r the scalar field $. The parameters 
are set to k = 1, z — 3, A = 23, u = 3.1215. 



As for massive bound KK modes, they are mainly impacted by two important param- 
eters, i.e., the coupling constant A and the position of the thin brane zq. For simplicity, we 
set it = in the following discussion. 

From the shapes in Fig. |9|, it is clear that the depth of both the potential well and 
the potential barrier increases with the coupling constant A, which means that there may 
exist more massive bound KK modes and resonances for stronger coupling. For example, 
for k = 1,Zq = 3 and u = 0, the numerical calculation shows that, there exists only one 
massive bound KK mode when A = 23 or A = 43, and exist two massive bound KK modes 
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when A = 53. The mass spectra are listed as follows: 



m 2 n = {9.31} U [21.53,oo) 
ml = {13.54} U [40.25,oo) 



71 



in 



n 



{15.29,46.27} U [49.61, 00) 



for A = 23, 
for A = 43, 
for A = 53. 



(3.41) 



And when A = 23, there is no resonance, but when A = 43, there are three resonances 
with m 2 = 41.0101, m 2 = 62.96 and m 2 = 78.2. The shapes of them are shown in Fig. 11. 
When A = 53, we found two resonances with m 2 = 71.829 and m 2 = 91.1831, the shapes 
of their wave function are shown in Fig. 
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Figure 11: (Color online) The shapes of the scalar resonances. The parameters are set to k = 
1, zq = 3, u = 0, A = 43. 

As shown in Fig. |9|c, with the increase of the distance of the two thin branes, the value 
of the potential Vq{zq) at z > zq will first increase and then decrease, which also can be 
seen from Fig. where we plot the shape of Vq{zq) as a function of Zq. The maximum of 
Vq{zq) can be calculated from Eq. (|3.40| ). The zq corresponding to the maximum of Vo(zq) 
is the positive root of the equation d ZQ Vo(zo) = 0, i.e., 



u 2 kzo = 6A arctan(fezo) (/czq arctan(fe2:o) — l) 



(3.42) 
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Figure 12: (Color online) The shapes of the scalar resonances. The parameters are set to k = 
1,zq = 3, u — 0, A = 53. 

Vo(zo) 




Figure 13: (Color online) The shape of the scalar potential Vq(zq) (as the function of zq). The 
parameters are set to k = 1, A = 23. 

Especially, for the case of u = 0, the above equation is reduced to l/(kzo) = arctan(/czo). 
The only one solution for this case is zq = 1. 16234 -r- This result shows that, when the 
distance of the two thin branes is of about the thickness of the thick brane, the potential 
that the scalar KK modes feel would be strongest, and hence there will be the most number 
of the massive bound KK modes. For example, the number of the massive bound KK modes 
with k = 1, Zq = 1.16234 is more than the case of k = 1, zq = 3 for A = 23, or A = 43, or 
A = 53, which can be seen from the following mass spectra: 

m 2 n = {9.31, 28.24, 41.12} U [43.45, oo), for z = 1.16, A = 23, 

ml = {13.56,41.01,63.00,78.30} U [81.23, oo), for z = 1.16, A = 43, (3.43) 

ml = {15.30, 46.27, 71.83, 91.23} U [100.12, oo), for z = 1.16, A = 53. 

However, when Vq{zq) reaches the maximum, there will be no resonances. What interesting 
is that the eigenvalues of some of these bound KK modes are almost that of the resonances 
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with some larger zq. For example, when k = 1, zq = 3, X = 43, there are three resonances 
with the eigenvalues m 2 = 41.0101, m 2 = 62.96 and m 2 = 78.2, and when k = 1,Zq = 
3, A = 53, there are two resonances with the eigenvalues m 2 = 71.829 and m 2 = 91.1831. It 
means that, when the distance of the two thin branes increases, some scalar bound modes 
would disappear, and the resonances are the remnants of the bound modes. 

3.3 Spin 1 vector fields 

Lastly, we will investigate the localization of spin 1 vectors. The action of a 5-dimensional 
vector field is 

Sl = ~\j d5x ^ GMNGRSF MRFNS, (3.44) 

where Fmn = 9m An — 9nAm- Then using the background metric (p^) , the equations of 
motion of this system are 

wv A FpA ) + rj^e- A d z (e A F 4A ) = 0, (3.45) 

rTdpFvi = 0. (3.46) 

Next considering the gauge invariant of § Ay[dx , we choose A4 = using gauge 
freedom. And with this gauge we make a decomposition of the vector field A a {x,z) = 
Yin a t l \ x )Pn( z ) e ~ A ^ 2 1 an d we can obtain that the KK modes of the vector field satisfy the 
following Schrodinger-like equation 

l-d 2 z + V^z)] p n (z) = m 2 nPn {z), (3.47) 

where m n are the masses of the 4-dimensional vectors, and V\(z) = 4[ — \- 4j-. Then 
substitute the Eq. ( 3.47 ) into the action ( |3.44 ), providing the orthonormality condition 



dz p m (z)p n (z) = 6 mn , (3-48) 
we obtain the 4-dimensional effective action: 

si = J2j d ' x ^a ( - i^r^/^/S - \< tr^f^) , (3.49) 

where g^ u = rf v , and f$ = d^a^ — d v affl is the 4-dimensional field strength tensor. 
For the GRS-inspired brane the effective potential of vector KK modes is read as 

TWl4(TTF^) -2;. W < z o {3 5Q) 

0, \z\ > Zq 
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The shapes of the potentials are shown below in Fig. 14. We can obtain the zero mode by 
solving Eq. ( 3.47| ) with m 2 = 0: 



p (z) oc 



[l + ikzfY 1 ' 4 , 

;i + (^ ) 2 )" 1/4 



\z\ < z 

\z\ > Zq 



(3.51) 



For this vector zero mode, the normalization condition ( 3.48| ) can not be satisfied both 
with a finite zq or an infinite zq, as the integral f p\dz always tends to infinity. So the 
massless vector can not be localized on the central thick brane. In fact, this result can be 
supported by the shape of the potential ( [3.50 ): it vanishes beyond the two thin branes for 
a finite zq, which indicates that there do not exist any bound KK modes including the zero 
mode. While for an infinite Zq, although the potential becomes a volcano potential, it is 
not enough to localize the zero mode. It is turned out that the result is the same as in the 
RS model case. It is known that in the RS model in AdSs space a spin 1 vector field is not 
localized neither on a brane with positive tension nor on a brane with negative tension, so 



the Dvali-Shifman mechanism should be considered for the vector field localization [59]. In 
Ref. |44] 1, it was also shown that a spin 1 vector can not be trapped on some Weyl thick 
branes with infinite fifth dimension obtained in Refs. [^, ^2], ^3|. However, it is turned 



out that it can be localized on the RS brane in some higher-dimensional cases [26|, or on 



a thick de Sitter brane [27 \ and on a Weyl thick brane with finite extra dimension 
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Figure 14: (Color online) The shapes of potentials of vectors V\. The parameters are set to 
k = 1, Zq = 3 for the thick line, and k = 1, Zq = 7 for the thin line. 



4. Discussions and conclusions 

In this paper, we have investigated the localization problem and mass spectra of various 
bulk matter fields such as spin ^ fermions, spin scalars and spin 1 vectors on the GRS- 



-22- 



inspired brane by presenting the mass-independent potentials of the KK modes in the 
corresponding Schrodinger equations. In the braneworld set-up, the extra dimension is 
infinite and non-compact, and there are one thick brane located at the origin of the extra 
dimension and two thin branes located at z = ±Zq. The bulk between the two thin branes 
is AdSs, while beyond the two thin branes it is Minkowski. Hence, the spectrum for all 
the fluctuations is made of continue KK modes (for volcano-like potentials) or discrete and 
continue KK modes (for Poschl- Teller- like potentials). 

For spin ^ fermions the formation of the potentials (3.6a) and ( |3.6b| ) has two sources: 
the gravity-fermion coupling $r M WM^ and the scalar-fermion coupling r]^>F((j))^ . If the 
space-time is flat (e A = 1), namely without the gravity-fermion coupling, the potentials will 
be Poschl- Teller-like ones for a kink solution, then the massless mode of left or right chiral 
fermion can be localized on the domain wall without any other condition. While without 
the scalar-fermion coupling, there will be no bound KK mode because the potentials vanish. 
So in order to localize 4-dimensional fermions, some kinds of scalar-fermion coupling should 
be introduced. 

In the GRS-inspired model, the potentials for both chiral fermion KK modes are pos- 
itive constant for the scalar-fermion coupling with F((ft) = 4P beyond the two thin branes. 
While between the two thin branes, where the bulk is AdSs, the minimum values of the 
potentials for left chiral and right chiral fermions are negative and positive for positive 
coupling constant 77, respectively. Hence, we always have massless left chiral fermions on 
the thick brane for a finite zq and a positive 77. But for an infinite zq, only when the 
coupling constant is larger that the critical one, the zero mode for the left chiral fermion 
can be localized. Furthermore, there exist some discrete bound fermion KK modes and 
continuous unbound ones. The number of the massive bound KK modes is impacted by 
three main parameters r],p,ZQ. When the scalar-fermion coupling constant r] or the cou- 
pling strength parameter p increases, the depth of the potential increases, which results in 
more 4-dimensional massive Dirac fermions trapped on the thick brane. It was also found 
that, when the distance of the two thin branes is not too small and the coupling constant 
T) is large, there exist resonance fermions. 

For spin scalars we have considered two cases, the free massless scalars and 
the massive scalar coupled with itself and the background scalar field <fi via a potential 
V($,4>) = (|^0 2 ~~ \ u2 ) ^ 2 ■ F° r t ne free massless scalars, the value of the potential Vq{z) 
is zero at \z\ > zq. So, the zero mode can not be localized on the thick brane if the position 
of the two thin branes zq is finite. However, if zq tends to infinity, the zero mode can be 
localized. When the distance of the two thin branes is much less than the thickness of the 
thick brane, there are resonances with spectrum m 2 oc n 2 , which is similar to the one of 
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the corresponding infinite square potential well. 

For the massive scalars, we need a fine tuning condition to ensure the localization of 
the scalar zero mode. With the fine tuning condition, one can make sure that there is no 
unstable scalar KK mode in the model. The number of the massive bound KK modes is 
mainly impacted by the parameters A and zq, which is similar to the case of fermions. With 
the increase of the coupling constant A, more massive bound KK modes would appear. We 
also found that there are scalar resonances under some conditions. 

For spin 1 vectors, there is no bound KK mode because the effective potential felt by 
vectors vanishes outside the two thin branes. 

In brief, in this braneworld model with the GRS geometry, the effective Schrodinger 
potentials for the fermions take the form of a modified Poschl- Teller potential, which in 
particular has a finite barrier, for the Yukawa and generalized Yukawa couplings to the 
background kink scalar. So there exist both bound KK modes and resonances, which are 
different with that for the potentials closer to the standard Poschl- Teller potentials without 
the finite barriers in Refs. [27, 29]. In fact, it can be seen from ( |3.6[ ) that the mass 
spectrum is decided by three factors: the warp factor e A , the scalar <fi, and the coupling 
F(4>), namely, the result in the paper is due to a combination of the GRS geometry, the 
kink like scalar, and the type of couplings to this background that have been chosen. For 
the scalars, with the coupling V(&,(fi), the effective Schrodinger potential is also modified 
Poschl- Teller potential with a finite barrier, which leads to the existence of both the bound 
KK modes and resonances. While in Refs. [ )56| | and [[58] with the AdS warped geometry, 
the Schrodinger potential for the scalars with the same coupling V(3>, (ft) is a volcano-like 
potential or an infinite potential. 
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